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Abstract 

We have calculated all the components of the current in a short one- 
dimensional channel between two superconductors for arbitrary voltages and 
transparencies D of the channel. We demonstrate that in the ballistic limit 
(D ~ 1), the crossover between the quasistationary evolution of the Joseph- 
son phase difference ip at small voltages and transport by multiple Andreev 
reflections at larger voltages can be described as the Landau-Zener transition 
induced by finite reflection in the channel. For perfect transmission and van- 
ishing energy relaxation rate the stationary current-phase relation is never 

recovered, and I(<p) = I c | simp/ 2 | signl/ for arbitrary small voltages. 
PACS numbers: 74.50.+r, 74.80.Fp, 73.20 Dx 
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It has been known for more than 20 years that electron transport in short superconduct- 
ing weak links with arbitrary transparency can be described in terms of multiple Andreev 
reflections (MAR). |]J Despite this, quantitative understanding of the ac Josephson effect in 
these structures is still not complete. Various approaches to quantitative calculations of the 
current at finite voltages were mainly focused on the dc current which exhibits the so- 
called subharmonic structure, i.e. current singularities at voltages V = 2A/en, n = 1,2, 
where A is the superconducting energy gap - see, e.g., @,|7j and references therein. However, 
dc current carries only indirect information about weak link dynamics, whereas calculations 
of the ac currents have been limited to large voltages, the limitation being caused by 
the fact that at small voltages it is necessary to take into account increasingly large number 
of Andreev reflections. 

The aim of our work was to study a model of a short constriction between two super- 
conductors which permits the quantitative description of the current dynamics for arbitrary 
voltages and transparencies of the constriction. In this model, we found a new regime in the 
constriction dynamics, which occurs at small voltages and connects quasistatic variations of 
Josephson phase difference at V — > with MAR at larger voltages. 

We consider a single-mode channel of electron gas with transparency D between two 
superconductors (the calculations can be generalized in a straightforward way to several 
separable modes). The length d of the channel is assumed to be much smaller than the 
coherence length £ as well as elastic and inelastic scattering lengths in the superconductors. 
This allows us to neglect scattering in the vicinity of the channel (besides that described 
by the reflection probability R = 1 — D) and makes it convenient to describe electron 
motion in the constriction with the time- dependent Bogolyubov-de Gennes (BdG) equations. 
Assuming that the Fermi energy in the constriction is much larger than the energy gap A 
we simplify these equations further by adopting the quasiclassical approximation. Condition 
d <C £ makes the superconducting properties of the constriction itself irrelevant (even if there 
is a finite A in the constriction we can neglect it in the BdG equations on the small space 
scale given by d) (TIJ. It is easier to visualize electron motion in the channel assuming that 
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the constriction is normal (A = 0), so that the transport through the resulting SNS structure 
can be described directly in terms of the Andreev reflection at the two NS interfaces. We 
adhere to this framework in what follows. 

The final simplification is that impedance of the single- (or few-) mode channel is on the 
order of h/e 2 and is much larger than characteristic impedance of a typical external circuit. 
This eliminates the necessity (essential for a realistic description of the Josephson junctions 
with low resistance) of determining the dynamics of the Josephson phase difference and 
voltage across the channel self-consistently. We assume that the voltage is constant in time. 



The model we obtain is directly applicable to the atomic-size Josephson junctions fll 
which exhibit ballistic quantization of the stationary critical current ||12|| . Another context 
of current interest in which the model is relevant is high-critical-current Josephson junc- 
tions @,|13[ which are believed to be adequately represented as an ensemble of atomic-size 
microconstrictions each of which carries a few conducting electron modes. 

Under the assumptions outlined above, the BdG equations for transport in the constric- 
tion can be solved in terms of the two scattering processes for electrons and holes, along the 
same lines as in the stationary case |TJ]]. One process is Andreev reflection at the NS inter- 



faces characterized by the reflection amplitude a as a function of the quasiparticle energy 



e: 



1 e-sign^-A 2 )^, |e|> A, 
[ e-^-e 2 ) 1 / 2 , | e |< A . 

Another process is electron scattering in the constriction characterized by a scattering ma- 
trix: 



(r t \ 
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where | t | 2 = D and | r | 2 = R. The scattering matrix for holes is the time reverse of S e i, 
Sh = S*i. 

The last ingredient of the scattering scheme is the fact that the energy of an electron is 
increased by eV each time it passes through the channel from left to right, while the hole 



increases its energy passing through the constriction in the opposite direction. Because of 
this the electron and hole wave functions are sums of the components with different energies 
shifted by 2eV. For instance, the wave functions in region I (Fig. 1) generated by the 
quasiparticle incident from the left superconductor onto the channel can be written as: 

i>A = E[( fl 2nA„ + J6 n0 )e ikx + Bn e- ikx ]e-^ +2neV ^ h , 



ij h = Y,[A n e lkx + a 2n B n e- ikx }e 



ikxi, —i(e+2neV)t/h 



(3) 



where k and e are momentum (equal to the Fermi momentum) and energy of the incident 
quasiparticle, and a m = a(e + meV). In eq. @ we took into account the fact that the 
amplitudes of electron and hole waves are related by the Andreev reflection, and the quasi- 
particle incident from the superconductor produces an electron in the normal region with 
amplitude J(e) = (1— | a(e) I 2 ) 1 / 2 . The wave function in region II has a similar form with 
two modifications: it does not have the source term J and is shifted in energy by eV. 
The wave amplitudes in regions I and II are related by the scattering matrix (0): 
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Eliminating the wave amplitudes C n , D n in region II from eq. @ we obtain the recurrence 
relation for the amplitudes A n , B n : 
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These recurrence relations can be solved with the method developed in [13,3]. The 
amplitudes A n , B n of the wave functions (^) obtained in this way determine all Fourier 
components of the current I(t) in the channel: 



I(t)=J2l k e 

k 



i2keVt/h 



Collecting contributions from the quasiparticles incident on the channel from the two su- 
perconductors, and making use of the fact that A(— e, —V) = — A*(e, V) and B(—e, — V) = 



Some results of the numerical calculations of the current from the recurrence relations @ 
and equations (§) are shown in Fig. 2. One can see that ac components of the current exhibit 
the subgap singularities at V = 2A/n similar to those in the dc current. It is clear from Fig. 



the tunnel Hamiltonian expressions for the current. In the case of zero temperature (shown in 
Fig. 2) the only component of the current that is non-vanishing at finite reflection probability 
and small voltages is the stationary Josephson current. In particular, it can be checked that 
the limiting (V = 0) values of the sine component (Fig. 2b) coincide with the first Fourier 
harmonics of the stationary Josephson current I ((f) = (eDA/2h) sin<£>/(l — D sin 2 (y9/2)) 1 / 2 . 

Figure 2c shows that for large transmission probabilities D cosine component of the cur- 
rent is negative in the wide range of voltages below the gap voltage 2A/e. (Calculations 
based on eqs. at finite temperatures show that this feature is preserved at temper- 

atures up to about 0.3A.) This fact provides a possible resolution of the long-standing 
problem of sign of the cosine component of the ac Josephson current in tunnel junctions 
(see, e.g., |16|), if one assumes that due to non-uniformity of tunnel barriers realistic tunnel 
junctions always contain regions with high transparency. 

The feature of the curves in Fig. 2 which to our knowledge has never before been dis- 
cussed is the rapid variation of all current components at small voltages and small reflection 
probabilities. In order to understand this new feature we consider first the case of perfect 
transmission, D — 1. In this case the recurrence relations can be solved explicitly, 




2 (and straightforward to show analytically) that in the limit D — > eqs. (|5]), (||) reproduce 
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2n-l 

B n = 0, A n = J Y[ clj . (7) 

3=1 

This solution shows that in the limit V — > when electrons and holes increase their energy in 
Andreev reflections cycles very slowly, the amplitudes A n decrease very rapidly as functions 
of energy outside the energy gap, since | a(e) |< 1 at these energies. At the same time, the 
number of Andreev reflections inside the gap increases. This means that, in this regime, the 
dominant contribution to the current comes from the gap interval | e |< A, where A n are 
not decaying (a(e) = exp{— i arccos(e/A)}, | a |= 1) and the number of Andreev reflection 
cycles diverges as A/eV. This divergence is regularized by the fact that the particles getting 
into the energy gap originate only from a small energy range near the gap edges. With this 
understanding eq. (||) for the kth Fourier component of the current can be simplified for 
y C A/e as follows: 

eA A r 1 

Ik = — — tanh(— ) J dzexp{2ik arccosz}. (8) 

Equation (|8]) means that the current-phase relation at finite voltages is: 

I(<p) = I c | sin cp/2 | signV , J c = ^- tanh( — ) , tp = 2eVt/h , (9) 

in particular, the dc current is 2I c /ir, the sine part of the first Fourier component of the ac 
current is vanishing, while the cosine part is — 2/ c /37r. 

Expression (^) has a natural interpretation in terms of the quasistationary discrete states 
inside the gap that are responsible for the stationary Josephson current [1^,0]. Two such 
states with energies E± = Acosyj/2 (Fig. 3) carry, respectively, forward and backward cur- 
rent, and in the stationary case (ip = const) are occupied according to equilibrium Boltzman 
distribution. At finite voltage the energy of these states is changing in time due to evolution 
of ip. For vanishing inelastic scattering rate the density of states in the superconductors 
is also vanishing within the gap, so that the occupation probabilities of the two current- 
carrying states remain constant as long as they are moving inside the energy gap, E± < A. 
The only point at which the occupation probabilities can change is at the gap edges e = ±A, 
(i.e. at <p = 27m). These considerations immediately give eq. (fj|). 



Solution (0) of the scattering problem at D = 1 gives some analytical results at larger 
voltages also. In particular it can be shown from eqs. ©,(0) that at T = the sine part of 
the first Fourier component of the ac Josephson current vanishes identically at all voltages, 
while cosine part approaches the asymptote — 7 c (Aln2/2eV) at eV ^> A. Both of these 
results agree with the numerical results shown in Fig. 2. 

The reasoning that lead to eq. @ can be easily generalized to the small but finite 
reflection probability R of the constriction. Finite reflection creates a finite matrix element 
r of the transition between two current- carrying states in the energy gap which occur near 
the point tp — 7T, where the energies E± of these states coincide -see Fig. 3. The problem 
of this transition is then a standard level-crossing problem and the probability p that the 
system will continue to occupy the same level after crossing the point ip = 7r is given by the 
Landau-Zener expression. In our notations this expression is: 

r tt-RA . , 

p = exp{ — ~^-\- ( 10 ) 

The finite transition probability p modifies the current-phase relation (|]) as follows: 

{sin 09/2 , < ip < 7T , 

(11) 
(2p - 1) sin<y?/2 , 7c<Lp<2n. 

In the relevant range of parameters (small R and V) eq. ( |TT1) reproduces the result of 
numerical solution of the recurrence relations (||]). Indeed, at very small V (eV < RA), p — > 
and the system in its evolution follows a current-phase relation which at T = coincides with 
the stationary relation I = I c sin p/2 sign(cos (p/2). As a result, the dc current and the cosine 
component of the first harmonic of the ac current are vanishing, while sine component is 
equal to its stationary value, in agreement with Fig. 2. At larger voltages (eV > -RA), p — > 1 
and the current-phase relation approaches the one for D = 1, 1(f) = JcSignV^ | sinip/2 |. 
For this 1(f) the sine component is zero while cosine component and dc current are non- 
vanishing. All this means that the reason for the rapid variation of all current components 
with voltage at D ~ 1 and small voltages (see Fig. 2) is that the probability (|l0| ) of Landau- 
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Zener transition between the two current-carrying states changes rapidly on a small voltage 
scale given by RA. 

Before concluding, we would like to mention that our calculations agree with most of the 
previous results on the ac Josephson effect in short constrictions, in the parameter ranges 
where previous results are available. In particular, our numerical results (Fig. 2) agree with 
the numerical results of Arnold || for large voltages. At D = 1 eq. (|9]) of our work gives 
the dc current in agreement with that obtained by Gunsenheimer and Zaikin ||. There is, 
however, a contradiction between calculations at large voltages based on the solution (0) 
and Zaitzev's results for ac components of the current at D — 1 and large voltages |J. The 
reason for this contradiction is not clear at present. 

In conclusion, we have calculated the current in a short single-mode electron channel 
between two superconductors for arbitrary voltages and transparencies of the channel. To 
the best of our knowledge this is the first time when full description of the current dynamics 
in a weak link is developed. In the ballistic limit, D ~ 1 crossover from quasistationary 
Josephson current at smaller voltages to multiple Andreev reflections at larger voltages 
occurs at V ~ irRA/e and can be described in terms of Landau-Zener tunneling between 
the discrete current-carrying states in the energy gap which are responsible for the stationary 
Josephson current at V = 0. 

The authors gratefully acknowledge numerous useful discussions with M. Kupriyanov and 
K. Likharev. This work was supported by DOD URI through AFOSR Grant # F49620-92- 
J0508. 
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FIGURES 



Figure 1. Schematic energy diagram of a short one-dimensional channel with arbitrary 
transparency between two superconductors. I and II denote the portions of the channel 
separated by the scattering region. 

Figure 2. DC current, and sine and cosine parts of the first Fourier component of the ac 
current in a short quantum channel between the two superconductors at zero temperature. 
Gt is the normal-state conductance of the channel, Gt = e 2 D/it%. All curves exhibit 
subharmonic singularities at V = 2A/n associated with multiple Andreev reflections. Rapid 
variation of the curves with D ~ 1 at small voltages is a manifestation of the Landau-Zener 
transitions between the quasistationary current-carrying states in the energy gap. For details 
see text. 

Figure 3. The energy diagram of the two quasistationary current-carrying states E± 
in the constriction. The arrows show two possible routes of the system evolution due to 
Landau-Zener transitions in the vicinity of the level-crossing point ip = ir. 
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